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Constraining the New Aether: Gravitational 
Cherenkov Radiation 



Abstract: We study the simplest concrete theory for spontaneous Lorentz viola- 
tion, the "New Aether Theory" of Jacobson and Mattingly, which is a vector-tensor 
gravitational theory with a fixed-modulus condition on the vector field. We show 
that the observation of ultra-high energy cosmic rays (which implies the absence of 
energy loss via various Cherenkov type processes) places constraints on the param- 
eters of this theory, which are much stronger than those previously found in the 
literature and are also stronger than the constraints generically arising when gravity 
displays sub-luminal propagation. 
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1. Introduction 

Lorentz invariance is a cornerstone of modern particle physics and of general relativ- 
ity. Exactly how well constrained is violation of Lorentz invariance? This question 
has received increasing experimental [Q and theoretical scrutiny recently (for a 
nice review see 0). Partly this is because one should always test the foundations of 
theories. Partly it is also because Lorentz violation is one possible route by which 
very high energy (Planck scale) physics might be detectable at the energy scales ac- 
tually available experimentally; we can constrain (rule out or perhaps find evidence 
for) theories involving the highest energy scales if those theories happen to predict 
any deviation from Lorentz invariance. The study of Lorentz invariance violation 
has also been stimulated by the apparent non-observation of the GZK cutoff in the 
cosmic ray flux (for a recent review see |^). 

Constraints on Lorentz violating physics are already quite severe. In fact, if 
Lorentz invariance is not a symmetry at the fundamental scale, then it is generi- 
cally communicated to low energy scales in an unsuppressed manner P, 0, leading 
to O(10~^) violations in experiments which instead provide < 10~^^ bounds on its 
violation. Similarly, constraints on high dimension Lorentz violating operators show 
that they are excluded if they arise at the Planck scale, even with coeflicients sup- 
pressed by several orders of magnitude M, (Supersymmetry provides a potential 



way around this problem |Tn . 
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Therefore, the study of Lorentz violation has shifted towards investigating the 
spontaneous violation of Lorentz invariance, arising at a scale smaller than the Planck 
scale. If this violation occurs in a way which distinguishes between different species 
of matter, then the constraints again become very severe; so one line of work has 
concentrated on spontaneous Lorentz violation due to a field coupled only to gravi- 
tation. 

Spontaneous Lorentz violation will occur whenever a field, transforming nontriv- 
ially under the Lorentz group, takes on a vacuum expectation value. This is the 
case, for instance, in the vector-tensor theories originally investigated by Will and 
Nordtvedt in the early 1970's Recently, a new twist on this idea, called the 



"New Aether Theory," has been proposed (|T^, [T3|, see also |TJ]). In this theory, 
besides the metric g^i, and ordinary matter fields, there is a new vector field S'^, 
which is constrained, either by a potential term or by a Lagrange multiplier, to take 
on a fixed length. The most general Lagrangian up to second order in derivatives for 
such a theory is, 

S = I d^x^g {-R -Cv- KgapS'^S^ - + Cm) , (1.1) 

IbTrGN J 

with Cm the matter Lagrangian (for instance, of the Standard Model) and 
Cv = ^ V„^^V"^^ + ^V^S'^VfsS^ + ^V„^^V^^" + ^^"V^^'^^^V^^^ , (1.2) 

yZ yZ ^ yZ yi 

which is the most general kinetic term quadratic in derivatives. (We follow the 



notational conventions of [0, for the Cj, which is opposite that of |T5[.) The 
potential term will ensure that the vector field S will take a nonzero timelike vacuum 
expectation value, S"^ = (we use [-1-,—,—,—] metric conventions), which then 
selects a preferred frame, the frame where the expectation value is purely in the time 
direction. In the limit A oo, this is a strict constraint; at finite but large A, the 
"radial" fluctuation is heavy (though we will argue below that the finite A theory is 
not viable). 

This "New Aether" theory is hardly the unique theory for studying the sponta- 
neous violation of Lorentz symmetry; but it is the simplest, and its study may teach 
us about what physical phenomena best constrain such theories, and what physical 
problems are encountered in even trying to formulate them. Therefore we consider it 
an interesting problem to study what the tightest physical limits on this theory are. 
Previous studies have focused on cosmological tests |]l6l, with bounds of order 10~^ 



on the Cj, and on the contribution of 5* to the tightly constrained post-Newtonian 
parameter a2 which puts a constraint of order 10~^ on a complicated function of 



the Cj, constraining the Cj to be 0(10^'') (except in special parameter ranges, where 
there can be cancellations). However, we expect that a much stronger constraint 
can be placed on this theory by calculating how efficiently relativistic particles will 
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radiate gravitational and S-field excitations. One then uses the arrival of very high 
energy cosmic rays to place a bound on how efficient such energy loss can be, and 
therefore on the theory's parameters. Similar considerations provide the tightest con- 
straints on subluminal propagation of gravity [^. However, the application to this 
problem is richer because there are five propagating modes, corresponding roughly 
to the two transverse gravitational modes and two transverse and one longitudinal 
mode of the S field fl^, and any of these can be emitted by a high energy cosmic 
ray. 

The purpose of this paper is to consider such processes and to use them to set a 
limit on the coefficients of the New Aether theory. As we shall see, the constraints 
are extremely strong, and represent the most stringent restrictions on this theory to 
date. Presumably this is because, while previous studies of this theory ||12|, |13|, |T6|, |I5| 
have generally been based on considering it as a classical theory, our treatment takes 
the theory to be valid semi-classically (in particle physics language, at tree level). 
While Einsteinian gravity is not consistent as a fully quantum theory, it is consistent 
at this same semi-classical level, so we consider it reasonable to ask for the New 
Aether theory to be interpreted as a semi-classical, and not just classical, theory. 

The outline of the remainder of the paper is as follows. In Section H, we make 
some remarks about the model, including bounds arising from causality, and the dif- 
ficulties encountered above the scale where radial excitations occur. In Section ^ we 
derive the constraint on the coefficients arising from the emission of traceless-tensor 
modes from high energy hadrons; Section ^ obtains a similar constraint when the 
"longitudinal S field" mode is considered, and Section || obtains another constraint 
from the emission of pairs of such excitations-a constraint which persists in the limit 
that all Cj are taken to zero at fixed ratio. Finally, we present a summary of our 
findings in the discussion section. 



2. Remarks on the New Aether Theory 

Before discussing constraints on the New Aether theory arising from a specific pro- 
cess, we will briefly mention some structural constraints on the theory. 

To get a taste for the problems at hand, consider the theory in the absence of 
the potential term, A = 0, and without a vacuum expectation value (using q rather 
than Cj/u^ as coefficients on the different terms). In this case, the back- reaction on 
the metric can be neglected. Since we are also temporarily not interested in the case 
with nonzero vacuum value, we can neglect the C4 term in Eq. (|1.2|) . The resulting 
effective Lagrangian is, 

Ar = -c^d^S^d'-S" - ic2+cs)d^S,d''S>' , (2.1) 

where we replaced covariant with ordinary derivatives, since the metric is taken to 
be fiat, and integrated by parts to make the C2 term look like the C3 term. 
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The point is that, even at the classical level, this Lagrangian yields a theory with 
an energy which is not bounded from below, except in the special case C1+C2+C3 = 
0. Therefore, as soon as interactions of any sort are included, we expect runaway 
behavior. The fastest way to see this is to consider the case where C2 = C3 = so 
only the first term is present; then express it in non-covariant language as 

£ = -ci (^d^Sjd^Sj - doSjdoSj - diSod^So + doSodoSo) . (2.2) 

This looks like an ordinary kinetic term for three fields named Si, S2, and 5*3, and 
a wrong-sign kinetic term for a fourth field named Sq. The associated Hamiltonian 
has positive energies associated with gradients and time derivatives of the Si, but 
negative energies associated with Sq: 

H = ci (d,S,d,S, + doS.doSj - d,Sod,So - 5o5o5o5o) . (2.3) 

This is the Hamiltonian density derived by finding canonical momenta for each S"^ 
field and applying a Legendre transform, but it does not coincide with the Hamilto- 
nian density Tqq derived by variation with respect to g^^. However, the difference is 
a total derivative, the total energy remains unbounded from below, and the theory 
is unstable. 

The contribution to Tqo arising from the (C2+C3) term also leads to an unstable 
theory, but as is well known, the instabilities cancel and the Hamiltonian is positive 
semidefinite for the special "gauge invariant" case (C1+C2+C3) = 0. In fact, this is the 
reason that, when discussing vector fields in flat-space quantum field theory, gauge 
invariance is required in any renormalizable, interacting theory. 

The theory with a vacuum expectation value but without a potential is also un- 
stable except in the gauge invariant case. Classically, as soon as there are excitations 
in the gravitational and S field degrees of freedom, the nonlinear interactions will 
allow a runaway in which positive energy excitations combined with the negative en- 
ergy excitations associated with the wrong-sign component of S grow without bound. 
In the quantum theory, even in vacuum and at tree level, there will be spontaneous 
production of a pair of positive and a pair of negative energy excitations; the rate 
is unbounded and this process can be used to rule out the theory, in analogy to 
"phantom" scalar theories ||T8[| . 

Adding a potential term —X{Sf^S^ — v^)^ with finite A does not rescue the theory. 
Expanding S^ = vu^ + 6S, the potential introduces a mass term for 6S^, 

Cd~4Xv\So)\ (2.4) 

but since the kinetic term for Sq has the wrong sign, this actually makes the field 
tachyonic. For instance, if C2 = = C3, the full Lagrangian for 5*0 is 

{-4Xv^)S^ + d^SodiSo - doSodoSo uj^ = -4Xv^ + . (2.5) 
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Changing the sign of A does not help because it destabihzes the Si fields, which then 
have a negative quartic term. 
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The only way out is to replace the potential with a strict constraint, S^S^^ = v 
This removes one degree of freedom (perturbatively, Sq) from the model, leaving 
a positive definite Hamiltonian. For this reason, in this paper we will take the 
New Aether theory to have the condition S^^S^ = imposed as a strict constraint. 
Of course, this nonlinear sigma model is not renormalizable, which means it lacks 
a good interpretation on scales larger than v beyond a semi-classical (tree) level. 
However, the same is true of gravity and the scale mpi, so we will not take this as a 
fundamental objection to considering this theory (though the lower v must be, the 
more of an embarrassment it becomes that the theory cannot be considered at the 
loop level above this scale). 

Wave mode squared speed polarization (when k = z) 

transverse, traceless metric 1/(1— ci— C3) /ii2, ^11 = — ^22 

transverse Aether ~ TTu^^ ( P^ ) — [(ci+cs)/ s]Si, i = 1,2 

trace - ^ee Eq. (j;^) hpo = -2 6S0, - 2 6So/s'' 

Table 1: Wave mode speeds and polarizations in the gauge h^i = Vi^i = 0, from [p^]. 

Jacobson and Mattingly have derived the normal modes of propagation of the 
theory, their speeds of propagation, and their admixture between graviton and 5* field 
in a particular gauge [|l3l. Their results are reproduced, with a few approximations, 
in Table ^ For future reference, the complete expressions for the transverse and 
trace propagation speeds are, 

s2 = 2ci -cl + cj 

^tran™ 2(ci + Ci){l - Ci - C3) ' ^ ' 



and 

[Ci + C2 + C3)(2 - Ci - C4) 
(Ci + C4)[2(l + C2f - (Ci + C2 + C3)(l + Ci + 2C2 + C3)] ' 

while the complete expressions for ha in the trace mode are 



2 ^ yi^i -r 1^2 -r L.3; yz.- 1.1- 

trace / ^ _L ^.^2 N /1 o M ' \ ■') 



hn = h22 = -(ci + c,)SSo , /i33 = ^i±^^l^^l±^25^o • (2.8) 

C1 + C2 + C3 

We end this section by mentioning what restrictions can be placed on the coef- 
ficients of this theory by demanding stability, the absence of ghosts, and causality. 



This has been addressed previously (see for instance ||T9[, [jT5[); and we have fairly 
little to add, except to note how easily the constraints follow from these speeds of 
propagation of the normal modes. As shown in the table, there are three distinct 
propagating modes with different dispersion relations. In the limit where the q 
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are small (the relevant limit, since large values are excluded), these can be under- 
stood as the standard transverse traceless graviton mode, the transverse S field mode 
6S ■ k = 0, and the longitudinal S field mode, SS || k (in the frame where the 5* field 
condensate is in the time direction, which we will use throughout). 

Gradient instabilities appear when a squared speed of propagation is negative; 
from the table, this imposes the constraints, 

>0, £i±^>o. (2.9) 



C1+C4 C1+C4 

Further, the kinetic (doSi)"^ term is wrong sign and the S field excitations are "ghost- 
like" (and ruled out by the arguments of [0) unless 

ci+C4>0, (2.10) 

and therefore ci > and (C1+C2+C3) > separately. 

We impose in addition the constraint that propagation speeds not be super- 
luminal. After all, we are considering a metric theory, which should be causal in the 
usual sense. This constraint imposes, 

<1, <1, ci+C3<0. (2.11) 



C1-I-C4 C1+C4 

More precisely, the first two constraints should be derived from Eq. ( p.6| ) and Eq. 
( P^.Tp . Lim has argued [|T^ that these constraints are essential, since their violation 
could allow propagation around closed curves. However they cannot be considered 
to be on as firm ground as the other constraints. 

The constraints we find will be on the sizes of the q. can be re-scaled to 
force f ^ = 1, which has been done in the previous literature. In this case, one would 
say that we are constraining the kinetic terms of S to be small. Alternately, we can 
rescale 5*^ to make the Cj/f^ be order 1, say, ci/f^ = 1/2, so the kinetic term for 
transverse excitations is canonically normalized; then the vacuum expectation value 
= 2ci is small, and we would say we are finding constraints on the size of the 
vacuum expectation value of the field 5*. We prefer this language, since it makes 
it more clear that what is being constrained is the size of the spontaneous Lorentz 
violation. However, previous authors have generally chosen = 1, and we will 
continue to define the c,- to follow their notation. 



3. Constraints from Gravi-Cherenkov Radiation 

All constraints we impose on the New Aether theory will arise by considering the pos- 
sibility that high energy cosmic rays should undergo a spontaneous emission process 
producing one or more excitations of the coupled gravity-S* field modes. That is, we 
shall consider diagrams of the form shown in Fig. |I], in which the solid line is either 
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Figure 1: Left: gravitational Cherenkov radiation. Right: emission of a pair of S" particles 
via a virtual graviton. 



a fermion or a gauge boson, with a dispersion relation we shall take to be lightlike, 
and the double-wiggly line is a graviton. First we consider the former diagram, with 
a graviton in the final state. 

Normally this process is impossible, because there is no way to conserve both the 
energy and the momentum of the final state particles and still have them on-shell. 
This would be true if they both satisfied massive or light-like dispersion relations. 
However, we are investigating a theory where Lorentz symmetry is broken, and the 
dispersion relations are different. Since in particular the transverse traceless graviton 
modes will now have sub-luminal dispersion, with fcg = /cq/s^, < 1 the propagation 
speed in the preferred frame, these modes carry more momentum than energy, and it 
is kinematically allowed to emit one. The larger a particle's energy in the preferred 
frame (assumed to be the frame defined by the microwave background radiation^), 
the higher the emission rate; so the arrival of extremely high energy cosmic rays sets 
a limit on how far from light-like the gravity dispersion relation can be. 

The problem of calculating the rate for this process has previously been solved 
by Moore and Nelson |T^, but we will repeat the derivation here giving somewhat 



more detail, as well as making a more careful treatment of the partonic structure of 
the cosmic ray hadrons which turns out to slightly strengthen the constraint derived. 

Consider the emission of gravitons by a fermionic parton inside of a hadron. 
(Since the energy scale involved in the problem is large, a partonic discussion is 
necessary, see for instance [Q). In the gauge used by Jacobson and Mattingly 
the emitted excitation is purely a graviton. The vertex between fermion and graviton 
is (see for instance pO[] ) 



'-^y^^^ [7,(Pi.+P2.) + l.iPi,+P2,) - 2r/,.(fi+n)] , (3.1) 



^If the preferred frame differs from the frame of the microwave background, our bounds actually 
become stronger. This is because high energy cosmic rays have been observed arriving from many 
directions. In any frame at a large boost 7 3> 1 with respect to the microwave frame, some of these 
high energy cosmic rays were more energetic than in our frame; our results should be applied to 
these cosmic rays. 
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where Pi and P2 are the 4-momenta of the two external fermion states. (Our notation 
will be that capital letters are 4-momenta, while lower case letters are the magnitudes 
of their spatial parts.) Since the matter Lagrangian is unchanged by including the 
S field, this vertex is the same; what is different is the dispersion relation of the 
final state graviton, which has maximum speed s — 1/ = 1 — e, with 
e ~ {—ci—cz)/2 <^ 1. (The assumption e <IC 1 is self-consistent, since we are deriving 
a constraint on e which will show that it is order 10~^^.) 

Labeling the incoming momentum P and the outgoing graviton and fermion 
momenta K and P — K, we then have = = (P — K)^ and = —2ek'^, where 
k = \k\ is the magnitude of the graviton spatial momentum in the preferred frame. 
It follows that P-K ^-ek"^. 

The total emission rate of gravitons from this parton is 

^^hlw? ^^'^t^^-^)'] 2^5[s-^kl-\kn Y: \M'\ , (3.2) 

where 5^ is the squared matrix element, summed over final and averaged over 
initial spins and polarizations, and we have used the energy-momentum conserving 
delta function to perform the integral over the fermionic final state momentum. Since 
the dispersion relation in the on-shell delta function for the graviton is only simple 
in the preferred frame, this is the best frame to evaluate the integral. The delta 
functions can be used to reduce the K integration to two integrals, over k and 4> a 
trivial azimuthal angle; 



1_ f d^K 
2p 



J (2^ 27r5[(P-X)^] 2nS[{l+2e)kl-k'] = dk d<f> , (3.3) 



with K'^ — —2ek'^ and pk — p-k — ek{p — k). 

— * 

The matrix element is easiest to evaluate in the frame where k defines the z axis; 
in this frame, if we choose p to lie in the x — z plane, p^. — 2ep{p — k). We find it 
most convenient to compute separately the contribution of the /in — /i22 and /ii2 + /i2i 
polarizations, with the former giving 

E l-^l' = ^^(4P.)' Tr (/'-^)7-f 7^ = l&7:G^^p{p-k){Ap'-Apk+e) , (3.4) 
and the latter giving 

E l-^l' = ^^(4p.)^ TV {f-Wr^y = IGTTGNeMP - k)k' ■ (3-5) 
Summing these, and then doing the trivial integral, the integrated rate is. 



327rG'Ne^ 



r dkp{p-k){2p^ - 2pk + k^) = -G^e^p^ . (3.6) 

JQ 2 
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Similarly, a gauge boson interacts via the vertex [pO | 



WnG^ _ P2Cm^,p<7 + Di,u,pa{Pi, P2) + Gauge dependent) , (3.7) 



where the gauge dependent part vanishes when contracted on external polarizations, 
are the graviton polarization indices, p is the Lorentz index for the boson with 
momentum Pi, cr is the Lorentz index for the boson with momentum P2, and 

Df,u,pa{Pl, P2) = Vpi'PlaP2p - [rippPlvP2p + rii^pPl<7P2u " ?7p<7^1^Ai. + (/^ ^ '^)I3.9) 

The kinematics for graviton emission from a gauge boson are the same as for a 
fermion, but the matrix element differs; we find the matrix element squared, summed 
over physical final polarizations and averaged over physical initial polarizations, is 

^ \M^\ = levrGNe' [{2p{p-k)+ef + k\2p-kf] , (3.10) 

leading to an emission rate of 

r = ^GNeV. (3.11) 

5 

The process discussed would cause ultra- high energy cosmic rays to lose energy. 
If the energy loss timescale is shorter than the propagation time across the galaxy 
~ 10 Kpc (the shortest scale on which we believe cosmic rays must propagate- 
arguably the highest energy cosmic rays propagate tens of megaparsecs) , then the 
cosmic rays would not reach us; so the time scale for the highest energy cosmic 
rays to propagate without substantial energy loss must be longer than 10 Kpc ~ 
lO^^GeV"^. The transverse momentum "kick" a parton experiences in emitting a 
graviton, found above, is of order ~ p^e; this is also the scale on which parton 
distribution functions should be evaluated. We will see in a moment that this is 
about ~ (3000 GeV)^; such a transverse momentum kick is enough to cause a 
hadron to fragment, so a single graviton emission should be taken to destroy the 
hadron, and the relevant rate is the total rate of graviton emission. This must be 
determined by summing over the partons in the hadron. We saw that the emission 
rate scales as p^, so the relevant rate is 

V^G..VY.jyM^{fl, (3.12) 

The right value of to evaluate the parton distribution functions f{x) is deter- 
mined self-consistently from the value of e saturating the bound; taking the hadron 
responsible for the highest energy cosmic ray to be a proton, and using CTEQ parton 
distribution functions |2^, we find ^ j x'^f{x) x (3/2 or 12/5) 0.025. The result 



is dominated by the fermion, even though (3/2) < (12/5), because the weighted 
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PDF is about 6 times larger for a quark. With this in mind we will neglect gluon 
contributions in the following sections. 

Considering that the highest energy cosmic ray which has been detected arriving 
at the earth had an energy ofp~3xl0^^ GeV ||22|, and requiring that it was able 
to propagate 10 Kpc, we find a bound, 

^^^c^-c,^^^^^_,,^ (3.13) 



This bound is somewhat stronger than the one quoted in [0, because we took a 
graviton emission to fragment the hadron, whereas their treatment just summed up 
the energy losses due to such emissions. Note that the emission rate depends on the 
departure of the dispersion relation from luminal, e, as e^. This should be compared 
to ordinary Cherenkov radiation, where the energy loss rate goes as e. This is a 
spin effect. The polarization vector/tensor is transverse and only sensitive to the 
transverse component of the momentum, which is 0{e^^'^). For a spin-1 particle this 
appears at first order in the amplitude, for a spin-2 particle it appears at second 
order. This will be important in what follows, because the S'-field polarization states 
are lower spin, and will therefore be less suppressed in the e — > limit. 



4. Emission of S modes 

The New Aether theory has five physical propagating modes, rather than the two 
familiar from pure general relativity. This can be simply understood; three degrees of 
freedom (the four S*^ modes, minus the mode removed by the fixed length condition) 
have been added. At small Cj, a gauge can be found in which the propagating modes 
divide into the two transverse traceless graviton modes and three 5* field modes. 

The approximate propagation speeds of these three modes are easily read off 
from the Lagrangian, Eq. ( |1.2|) , by replacing covariant with ordinary derivatives, 
and by setting S^V^ do in the last term; 

Cv = Si (^-^{kl - kiki6ij) '^^^2 ^^ kikj + -^^o^ Sj . (4.1) 

Two degenerate modes correspond to Siki = 0, the transverse Aether modes, with 
dispersion relation (ci/v^ + C4/v'^)kQ = {ci/v^)k'^. The remaining mode has Si \\ ki] 
we shall call it a longitudinal mode. Its dispersion relation is (ci/f^ + C4/v^)kl = 
{[ci+C2+Cs]/v'^)k^. These reproduce the small c limits presented in Table |I|. 

If we consider the f -scaled q to be 0(1) coefficients and <^ m'^i, then it is 
natural to treat C4 ~ cf. This interpretation is natural if the fundamental scale is 
the Planck scale but the vacuum expectation value is smaller. For this reason we 
will take the transverse mode to propagate close to the speed of light. Under this 
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approximation, we need the full expression, Eq. ( p.6|) , for its speed of propagation 



2 ^ 2ci -cj + cl _ + 2ciC3 + - 2C4 ^ . _ , .42^ 

Vansverse 2(01+04) (1 - Ci - C3) " 2Ci " 



Here e' ~ q is a small parameter indicating the failure of the transverse mode to 
move at the speed of light. 

When gravity is included, each mode takes on a small admixture of graviton, 
and therefore couples to standard model particles and can be emitted as radiation 
in processes analogous to the one considered in the last section. To compute the 
emission rate, we must include the leading order /i^^ content of the excitation to find 
the matrix element, and we must correctly account for the particle's spectral weight 
(correctly normalize the external state) in writing down the external phase space. 

We begin with the transverse S excitation. The correct on-shell delta function 
accounting for the particle's spectral weight (correctly normalizing the external state) 
is 

27r5(2{ci+Ci)kl - 2ci\k\'^^ , (4.3) 

which introduces a 1/ci factor in the rate. The vertex with a fermion is Eq. ( |3.1j ) 
contracted with hi^ = h^i = (ci+C3)/s (see Table |l|) for one polarization and /123 = 
^32 = {ci+cs)/s for the other. Therefore introduces two positive factors of 

(C1+C3) and the result is suppressed by c relative to the result for a graviton. However, 
as we will see, the matrix element will contribute only one, not two, factors of 1 — s^, 
so the rate is just as large as for a spin-2 graviton. 

Using the vertex and polarization just presented, choosing coordinates such that 
only is nonzero, the matrix element is. 



M = ^v/l6^(^i+^3) -(^_^^ ^^^(2P - K), + 73(2P - K),] u{P) {AA) 

for one polarization and 1 — > 2 for the other polarization. The spin summed, squared 
matrix element is 

^\M\^ = 47re'GN(ci+C3)^ (32/ - 32^^ A: + A2p^e - lOpfc^ + A;^) , (4.5) 

where e' is the correction from lightlike dispersion for the transverse mode, defined 
above. The total production rate from a fermion is obtained by integrating over k 
and is 

r = gI!(ci+C3)^(cf + 2c3Ci + ci-2c4) 3 . (4.6) 

80 2cf ^ ^ ' 

We have not computed the analogous expression for gluons because the parton dis- 
tribution function of gluons leads to a smaller contribution. 
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Using the same integral over parton distributions as in the last section and 
applying this bound to the highest energy cosmic ray observed, we get a constraint, 



(Ci+C3)2(c2 + 2ciC3 + Cj - 2C4) ^ ^ ^ ^^.33 



2ci 

Since the quantity constrained is ~ c^, this constraint is of comparable strength to 
the previous one-but has an independent functional dependence on the q. 

Now we turn to the longitudinal 5* mode. The dispersion relation for this mode, 
~ (ci+C2+C3)/(ci+C4), is the only one which generically deviates by an 0(1) 
amount from = 1 even in the limit of small q (small vacuum value v for the S 
field). It will therefore lead to a strong constraint. But as we will see, the constraint 
is even strong when ~ 1, because there is no spin suppression in the emission of 
this particle in this limit. 

We begin with the phase space. In the small q limit and in a gauge where the 
gauge condition is based only on h^j_y, the energy-momentum conserving delta func- 
tion correctly accounting for the spectral weight (right external state normalization) 
of the longitudinal mode is 

2Ti5(2{ci+c^)kl - 2(ci+C2+C3)|^p) , (4.8) 

leading to a phase space integration of 

^ / ^-^j^^2rf(P|)2^i [2(ci+C2+C3)F - 2(ci+C4)iS] (2ir)V(P-P,-K) 

J {2TTy A\Ci+C2+C3\k L 



2p 



1 /-ifip 



327r|ci-FC2+C3|p2 



dk, (4.9) 



where = {C1+C2+C3) / {C1+C4) is the propagation speed. 

However, Jacobson and Mattingly's results for the mixing of the graviton 



and S field in the state polarization are in a gauge where diSi = 0. To use the 
polarization vectors found there, we must relate the amplitude of 6S0 in this gauge 
to the amphtude of 55*^ in the above gauge. The relation turns out to be 

(55'o(Jacobson's gauge) = ^^^^(above gauge) . (4-10) 

Therefore the matrix element squared will contain an extra factor of if the matrix 
element is computed normalizing the vertex based on 5*0, which means |ci+C2+C3| 
should be replaced with |ci+C4|. 

The results we have written for the dispersion relation and polarization tensor 
in Table ^ are insufficient to compute the matrix element because it vanishes at this 
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order; 



u{P - K) 



7o(2P-ir)o-^73(2P-ir)3 



s2 



u{P) oc u{P - K) [i^.K^'] u{P) = . 

(4.11) 

Therefore we must use the full expressions for the dispersion relations and the polar- 
ization tensor, Eq. ( p.7| ) and Eq. (|2.8|), which expanded to the desired order become 
(choosing k as the z axis) 

^00 = -s'^hs - (02+04-03)5*0 , hu = = -(ci+C4)S'o . (4.12) 

These lead to a matrix element squared of 

J2\M\^ = 47rGN(2p - ks)\Ap^ - Aspk - k^l - s^)){cs-Ciy , (4.13) 

and a total rate of 

r = AG„p3 (fiz£li! fa' + 24.^ + 35. + 20 _ 

15 C1+C4 i^ + sr 



with s = (ci+C2+C3)/(ci+C4). If we had not expanded in small q, then (03— 04)^ 
above would instead read (03(2+04)— 2c4+ci(ci+C3+C4))^/(l—ci—C3)^. 

The function of s above ranges from 20 to 5.3125 as s goes from to 1; for 
simplicity we will set it to 5, which will give a conservative bound. Since p^. can 
be as large as for this process if s is small, we also use the parton distribution 
functions evaluated at 10^° GeV, which give J x^f{x) = .0052 for quarks. Using the 
same constraint on F as before, we find 

^^r:^ < 1 X 10-^\ (4.15) 
IC1+C4I 

Curiously, the emission rate which leads to this constraint does not vanish in the 
limit s — >^ 1 that the emitted mode approaches lightlike propagation. This is because 
its emission, unlike transverse S and graviton emission, is not spin suppressed. It 
can also be understood as an example of non-decoupling of a longitudinal mode of a 
vector field when the mass becomes zero but gauge invariance is broken. However, 
when (1 — s) < 10~^^, then it becomes important that the proton dispersion is not 
precisely lightlike; so the above limit is conditional on (c2 + C3 — C4)/ci > 10~^^. 

The emission rate also has a nonzero limit as s goes to zero, (C1+C2+C3) 0, 
which is the "gauge invariant" limit for the kinetic term under which the theory 
remains free of tachyons if the strict fixed-modulus limit on the S field is lifted. 
However, if the strict fixed modulus limit is lifted and Af^ < (10^° GeV)^, then the 
treatment of the external propagating states must be reconsidered. Note that, as s — > 
0, the energy lost by emitting a longitudinal S mode goes to zero. Naively, this means 
that the emission process does not degrade the energy of the emitter, and is therefore 
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allowed. The emission process can still be ruled out, however, because it causes the 
breakup of the hadron, so no single final state particle carries as much energy as the 
initial hadron, and the highest energy cosmic rays would still be degraded in energy.^ 



5. ?/' ^ iIjSS 



The other process we will consider is the emission of two S field excitations via an 
off-shell graviton propagator, the diagram on the right in Fig. ^ This process is sup- 
pressed, with respect to the processes we have just discussed, by an additional power 
of so naively it would not be important. However, the processes we considered 
all vanish as the Cj are taken to zero; the graviton dispersion approaches the light- 
cone, preventing its emission, and the direct coupling of matter to S field excitations 
(via their mixing with the graviton) is suppressed and also vanishes linearly in c. In 
contrast, the coupling of the graviton to matter and to (canonically normalized) S 
fields remains finite as the Cj are taken to zero. In this limit, the transverse S field 
excitation becomes lightlike and the phase space for its production disappears; but 
the longitudinal S excitation remains sub-luminal with = {C1+C2+C3) / (ci) and its 
production persists. Therefore it is possible to set a limit on (ratios of) the which 
persists even as the q are taken to zero (with ratios held fixed). 

In studying this process, we will systematically expand in small Cj, which is 
justified by the results of the previous sections. We will also take C4 ^ ci, which 
we previously argued is reasonable. Were we not to make this approximation, then 
transverse S modes could also be pair produced, and we would find that this process 
implies that c^/ci is indeed small. Since the q are small, we can systematically 
neglect the mixing of graviton and Aether excitations. 

We begin by making a crude parametric estimate of the production rate. At 
generic values of s^, the squared speed of the longitudinal mode, the production rate 
should be of order 



purely on parametric grounds. Substituting p ~ 10^^ GeV, this gives F ~ 10^^^ GeV. 
This is to be compared with the propagation time of a typical cosmic ray, ~ 10'^ Kpc 
~ lO'^^/GeV. The mean lifetime of a cosmic ray before emitting an S pair is much 
shorter than the propagation time. Therefore all parameters are excluded which do 
not make 1 — <^ 1, and we are justified to expand in this quantity. 

^If the highest energy cosmic rays are elementary particles such as photons, rather than hadrons, 
all of our bounds are much stronger, since J dxx^f{x) is approximately 1 in this case. For the 
s — > limit, energy loss occurs whenever the S emission is accompanied by bremsstrahlung or pair 
production. The inelastic cross-section for photon primaries is large enough that the energy loss 
rate still exceeds that for hadrons. 



F 
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5.1 Kinematics 

We begin with a treatment of the kinematics of the reaction. The theory has a 
preferred frame, and it is awkward not to work in this frame. Since this is not the 
rest frame of the emitting particle, and since the usual Lorentz invariance tricks for 
treating phase space cannot be used, the treatment of the kinematics is somewhat 
more complicated than would usually be the case for an integrated 3-body decay 
rate. We choose the coordinate axes such that the momentum of the initial particle 
is along the z-axis. The four-momenta of the 3 final-state particles are written as ki 
and k2 for the two longitudinal S modes and ka for the mode with lightlike dispersion. 
All but 5 final state variables are performed by energy-momentum conserving delta 
functions. We choose to leave as integration variables the following 5 independent 
variables: 

fci^, /C3^, Ika^l , P, and 0, 

where f3 is the angle between ki^ and — ks^. The angle corresponds to an overall 
rotation around the z-axis and cannot be fixed by the conservation of momentum. At 
the same time, this angle does not appear anywhere in the matrix elements involved 
in the emission processes, so the integral over this independent variable is trivial. 

We need expressions for the other variables in the problem in terms of these 
5 independent variables. However, we do not need absolutely general expressions, 
because we can make two approximations. First, defining 

s = l-e, (5.2) 

we have e -C 1. (The estimate, Eq. ( p.l|) , shows that e ~ 1 is already excluded, 
so we may take e <^ 1 in trying to determine the limit on e.) This means that, at 
generic kiz and k^z, the k are all coUinear, and we may expand to quadratic order in 
the k^. Second, the graviton propagator will have a virtuality of (P — K^)"^. When 
ka^ is small, this is close to zero, and the matrix element is enhanced by this nearly 
on-shell intermediate state. For this reason, the dominant region of the integration 
is at small ks^, and we may expand in this quantity being smaller than ki^ and k2±. 
Under the coUinear approximation, energy conservation reads 

1 \ J 2 f ^ 1 \ ^l±^3± COS P 

Kh. + P-k,z-kJ \m. ^ P-hz-kJ P-k,z-k,. -2e(p-A:3.) = . 

(5.3) 

Further approximating that ks^ ^ ki^, we find, 

kl^ = kzik^2V2e = kl^ . (5.4) 

To make the notation in the remainder of the section more compact, we will 
define x = , y = z = x, y, and z are independent integration variables, 
with X in the range [0, 1], ?/ in the range [0, x], and z ranging up to of order y/e, but 
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as discussed, the rate will be dominated by the region where z <^ ^e, so the upper 
bound of this integration will not be important. The /3 integral runs over the range 
[0, 2-7r] for small and at leading order in small z the integrand will not depend on 
/3, so this integral may also be conducted trivially. 

With these simplifications, the phase space becomes 



1 1 

2 2^ 
1 

210^ 
V 



(fki (fk2 cPh 
{2'KfSkik2h 
dkidks 



{27r)%k,, + k2, + hMp-kl~kl-kl)5\\^^^+'k2±+'k^±) 



dp ki±dki± k^^dk^^5 



dx I dy I dz 



x{l — x) 



k^ 



1 1 
ki k2 



- 2e(p - ks 



(5.5) 



where the small approximation was made in passing from the first to second ex- 
pressions. The leading factor of 1/2 in the first expression is the final state symmetry 
factor. 



5.2 Matrix Element 

The vertex coupling the graviton to two 5* field excitations can be found by expanding 
Eq. (|1.2| ) to linear order in hf^,^. As discussed above, we neglect C4. Writing the vertex 
as h^,{Q)S^{K,)Sp{K2)V^-'"P, we find 



-Cl 



-C2 
-C3 



^a(p^u)p ^2Ki ■ K2 + Q^) + 27]''^ K^^Ki' + 2{Ki-K2Y>'r]''^^'' Q^^ 
2K'^{K2 + QY^t]"^'^ - K^Q^r]^"' + 2K^{Ki + g)^^!]'')" - K^Q^'r]'"' 



This is linear in the q and so the matrix element squared is quadratic in the q; 
however, there is a l/(2ci) factor associated with each S field external state, see Eq. 



( [4.3|) , so the rate is zero order in the q. Alternatively, if we canonically normalize 
the S field from the start, each term in Eq. { f).6\ j should be multiplied by l/(2ci). 

In calculating the matrix element we must be somewhat careful, since the final 
state phase space remains non-vanishing even as the graviton becomes lightlike. Since 
the matrix element will be singular at this point due to the graviton propagator, 
one must treat the kinematics of the very small k^^ region carefully, in particular 
including the nonzero proton mass in computing Q^, how far off-shell the graviton 
is. (Because of this near-singularity, the dominant \Q^\ will be < 171^^ the proton 
mass squared. The process is therefore not probing the structure of the proton, and 
we can consider it as a single fermionic object.) This gives corrections of the form 
rripj./p'^ to the virtuality of the graviton. 
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We again use the expression for the graviton-matter couphng from and take 
the graviton propagator in the usual Dedonder gauge. The longitudinal polarization 
vectors are adequately approximated by 

ei ~ e2 ~ £ , (5.7) 

good at leading order in e^/^. Therefore, 

ei ■ p = 62 ■ p ^ p, 

el ■ 4 = £2 ■ 4 -P(l -a;). (5.8) 

Because of the near-singularity in the graviton propagator, the cross-section is 
dominated by the kinematic regime ~ < epg ~ i^Ti ■ 7^2 ~ ^ ■ ^i- In this 
regime most of the terms in Eq. (|5.6| ) are subdominant; only the terms involving 
ciKi ■ K2, C2{KfK^ + K^K^), and c^iK^K^ + K^K^) contribute at order epl 
Furthermore, since we are specializing to the regime e <^ 1, we must have IC2+C3I -C 
Ci; otherwise the propagation speed is far from luminal, which we have already 
excluded. Therefore the C2 and C3 contributions to the matrix element cancel in the 
leading kinematic regime and we are left with only the Ci contribution. This rather 
substantially simplifies the calculation. The matrix element, including the final state 
normalization factor mentioned above, reads 



(167rG 



{K, ■ K^f {z ■p){z- h) {z-p + z-h)\ (5.9) 



In terms of the variables defined in the preceding section this is 

2_^\M\ =£:(167rGN)p ^ -2 , (5.10) 



which is to be substituted into Eq. ( ^.5|) . The z integral gives, 

zdz 



(5.11) 



illustrating that the integral is dominated by small z as claimed. The remaining x 
and y integrations are straightforward, and give a total rate of 

V = lll^. (5.12) 



2257rmpj. 



Unlike the processes in the previous two sections, in this process the transverse 
momentum "kick" to the proton is relatively small, so the proton typically does not 
break up, and we should compute the rate of energy loss rather than the frequency 
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of emission. This is done by adding a factor xp inside the integral, leading to an 
energy loss per time of 

dt 144007rm2^ ' ^ ' ' 

Assuming that the energy loss should not exceed the total energy of the highest 
energy cosmic ray observed over a 10 Kpc distance, and again using p = 3 x 10^^ 
GeV this leads to a bound on e of 



e ^ ^i^^^^ < 3 X 10-1^ (5.14) 

Cl 

We see that the approximation {m^^/p^) ^ e is self-consistent. 
6. Discussion and Conclusions 

Spontaneous Lorentz violation requires that a field, transforming nontrivially under 
the Lorentz group, take on a vacuum expectation value. The simplest available field 
is a vector field, and the most general theory where such a field takes on a vacuum 
value, which is stable and free of ghosts, is the New Aether theory, with a Lagrangian 
given in Eq. ( |1 . 1| ) , Eq. ( |1.2| ). This theory is ill behaved unless the restriction on the 
modulus Sfj^S'^ is imposed as a constraint, not through a potential. This renders the 
theory non-renormalizable and valid only semi-classically. 

The most notable feature of the theory at the semi-classical level is that there 
are 5 propagating modes — two graviton, two transverse 5, and one longitudinal 5" — 
which are massless but have sub-luminal propagation. As a result, high energy 
particles moving close to the speed of light can radiate them, in analogy with the 
Cherenkov process. This makes high energy particles lose energy-the higher energy 
the faster the energy loss. Since very high energy cosmic rays are known to travel 
astronomical distances, this places constraints on the theory. Namely, it cannot have 
parameters causing energy loss which would degrade those high energy cosmic rays. 
Since there are several modes with different dispersion, there are several constraints. 
In particular, we have found, 

-Cl - C3 < 1 X 10"^^ Eq. ( Prri (6.1) 
(^^+^^^'(^' + 'r^ + ^'-2c4)^^_^^^Q_3i Eq.(i3) (6.2) 



< 1 X 10-3° Eq. (1^ (6.3) 



|Cl-FC4 

C4 - C2 - C3 ^ „ , „_i9 



< 3 X 10"^^ Eq. (|5:T4|) (6.4) 



Cl 

-22 



The next-to-last bound is conditional on (c4 — C2 — C3)/ci > 10^ . In every case 
the opposite side of the expression is bounded by 0. This is because no dispersion 
relation can be super-luminal in flat space. 



18 



Combining these bounds with Eq. (|2.9|) -Eq. the Cj are constrained such 

that each must be \ci\ < 10~^^, except for a special case in which each dispersion 
relation in Table |l| is almost exactly luminal. Using the full expressions, Eq. ( |2.6| ) 
and Eq. ( |2.7] ) for the propagation speeds, this special case occurs when the Cj satisfy 

C3 = -Cl , C4 = , C2 = ^\ . (6.5) 

i — Zci 

We have checked that these same conditions cause 02 as determined in to vanish 
(assuming the result quoted there is leading order in the Cj). 

Physically, this means that Cherenkov processes rule out the New Aether theory 
unless either 

1. parameters are just right such that all three propagating modes (graviton and 
two S field) are light-like to extremely high precision, a 1-dimensional subspace 
of the 4-dimensional parameter space of the model, or 

2. canonically normalizing the field S, the vacuum expectation value v of S sat- 
isfies f < 3 X lO^^mpi ~ 5 X 10^^ GeV, or 



the constraints, Eq. (|2.11|) , are not correct; there are modes which propagate 



super-luminally, despite the theory being a generally covariant metric theory 
of gravity. (We earlier discounted this possibility because we are considering 
a metric theory where the violation of Lorentz invariance is supposed to be 
spontaneous. However, the necessity of the constraint against super-luminal 
propagation is controversial, so we acknowledge this logical possibility.) 

The constraints are not fatal to the theory, but they mean that either the theory 
does not predict Lorentz violating propagation speeds, or the scale of the spontaneous 
Lorentz violation must be much smaller than the Planck scale. 
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